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Abstract 

This letter provides the first prediction of the probabihty density func- 
tion (PDF) of flux R in plasma sheath region in magnetic fusion de- 
vices which is characterized by dynamical equations with exponential 
non-linearities. By using a non-perturbative statistical theory (in- 
stantons), the PDF tails of first moment are shown to be a modified 
Gumbel distribution which represents a frequency distribution of the 
extreme values of the ensemble. The non-Gaussian PDF tails that may 
be enhanced over Gaussian predictions are the result of intermittency 
caused by short lived coherent structures (instantons). 
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In the edge region of magnetically confined fusion experiment the plasma 
exhibits large bursty fluctuations, often called intermittency [Tj- [9]. Inter- 
mittent turbulence is characterized by patchy spatial structure that is bursty 
in time. The likelihood of these events are described by probability den- 
sity functions (PDFs) which significantly deviate from Gaussian predictions. 
More specifically, exponential scalings are often observed in the PDF tails 
in magnetic confinement experiments [6]- [9], and intermittency at the edge 
strongly influences the overall global particle and heat transport. In par- 
ticular it may for instance influence the threshold for the high confinement 
mode (H-mode) in tokamak experiments [IQl. In the view of these experi- 
mental results, theories built on average transport coefficients and Gaussian 
statistics fall short in predicting vital transport processes. The high possibil- 
ity of confinement degradation by intermittency strongly calls for a powerful 
predictive theory. 

The purpose of this letter is to provide the theoretical prediction for the 
PDF tails of flux in plasma sheath by using a non-perturbative statistical 
method. The basic idea is to associate the bursty event with the creation 
of a coherent structure (e.g. blobs, streamers etc.). The candidate that 
could describe the creation process is the instanton which is localized in 
time, existing during the formation of the coherent structure. This idea is 
embedded in our theoretical method - the so-called instanton method - which 
is a non-perturbative way of calculating the PDF tails [11]- [22]. Note that, 
the instanton method has been adopted from quantum field theory and then 
modified to classical statistical physics for Burgers turbulence [12j- [13j and 
in the model by Kraichnan [23j. 

Our previous papers using this instanton method have shown that the 
PDF tails of momentum flux and heat flux are significantly enhanced over 
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the Gaussian prediction [T3]- [IB], providing a novel explanation for expo- 
nential PDF tails of momentum flux found in recent experiments at CSDX 
at UCSD [9j. PDF tails of forced shear flow i? in a non-linear diffusion model 
were shown to be ~ exp{— c(i?/-Ro)^} that falls off much faster than a Gaus- 
sian by both analytical (e.g. instanton) and numerical investigations [2T]. 

In the following we study the PDFs of flux of 1st moment variables such 
as potential, density etc in plasma sheath where electric potential is governed 
by exponential non-linear interaction. In Ref. [SU], we computed the PDF 
tails of blob density by using a simplified model for blobby transport with 
different non-linear interaction. Note that in the instanton analysis a careful 
treatment of the dominant non-linear interaction is necessary so that keeping 
only leading order non-linear interaction by Taylor expansion may lead to the 
erroneous predictions for the PDF tails [22j. By using the instanton method 
we show that the PDF tails of the first moment satisfy modified Gumbel 
distributions [23]- [25]. The Gumbel distribution is used to model extreme 
value distributions, i.e. a frequency distribution of the largest or smallest 
values of the ensemble. This is due to the fact that Taylor expansion of 
the exponential contains infinite sum of all powers. In a comparison with a 
Gaussian distribution it is shown that the flux in plasma sheath may exhibit 
stretched PDF tails. We emphasise that the predicted Gumbel distribution 
is new and has not been found in any related previous work. 

As a physical model of plasma sheath, we consider a generic simplified 
model of turbulence with a dominant exponential non-linear interaction. The 
model is a simplified version of the three field Braghinskii fluid model pre- 
sented in Ref. [26] where we have retained the exponential non-linear inter- 
action. Note that this model can be reduced to the model derived in Ref. [27] 
by linearization and keeping only non-linear terms up to second order. This 
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model is sufficient for the prediction of the PDF tails since our previous pa- 
pers showed that the PDF tails are rather insensitive to the details of the dy- 
namics and dependent only on the dominant non-linear interaction term [22] . 
We model the sheath region with a general non-linear equation describing the 
time evolution of the variable for fluctuations with exponential interaction 
Ke""^ 0, [271- [211, 

f + K^'' = /• (1) 

Here / is the forcing. The statistics of the forcing is assumed to be Gaussian 
with a short correlation time modeled by the delta function as 



and (/) = 0. The angular brackets denote the average over the statistics of 
the forcing /. The delta correlation in time is chosen for the sake of simplicity 
of the analysis. Note that in the case of a finite correlation time non-local 
integral equations in time are needed. 

We calculate the PDF tails of local flux M{(j){x = xq)) that is the ffist mo- 
ment of (f) (e.g. density, potential) by using the instanton method. The PDF 
tails are expressed in terms of a path-integral using the Gaussian statistics 
of the forcing [HI. The optimum path is then associated with the creation of 
a short lived coherent structure (among all possible paths) and the action is 
evaluated using the saddle-point method on the effective action. The saddle- 
point solution of the dynamical variable 0(x, t) of the form 0(x, t) = F{t)v{x) 
is called an instanton if F{t) = at t = — oo and F(t) 7^ at t = 0. Note 
that, the function v{x) here represents the spatial form of the coherent struc- 
ture. The probability density function of the flux R can be defined as 



{fix,t)fix\t')) = 6it-t')Kix-x'), 



(2) 




(3) 
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where 

7A = (exp(-iAM(0(x = xo)))). (4) 
The integrand can then be rewritten in the form of a path-integral as 

= J V(i)V4)e-^\ (5) 
Here the effective action S\ in Eq. (5) is expressed as, 
Sx ^ -ij dxdt^ + Xe^"^ j 

+ -J dxdx'dt4>{x, t)K{x — x')4>{x', t) 

+ i\ j dxdtM{(p)5{x - XQ)5{t). (6) 

We make use of a logarithmic transform = Inrt and u = v{F + 1) where v is 
the spatial structure of the instanton and F is the time dependent function. 
The action can now be recast into, 

= -i/*c.A(^M£±i)+c.(F+ir) 

+ ^ciC4 j dtFl 

+ ixJdtciCsln{F + l)S{t). (7) 
Here we have used the definitions. 



Ci = J dxv{x), (8) 

C1C2 = j dxv{x)Kv{xY, (9) 

C1C3 = j dxd{x — xq), (10) 

C1C4 = j dxdyv{x)v{y)K{x - y), (11) 

In order to find the path with highest probability identified by the instan- 
ton or the extremum of the action we compute the first variational derivatives 
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to minimize Sx with respect to F and Fi, 

^ = -i(^ln(F + l) + C2(F + ir)+C4A=0, (12) 

^ = -i(^-^ + C2c{F + ir'F,^+iXcsj^S{^ (13) 
Eq. (12) gives 

+ i [f + C2{F + iy+^) ^ C4{F + l)Fi (14) 
Differentiating Eq. (14) w.r.t time gives, 



F + C2(c + 1)(F + lyP - - C2(F + lyP^ 

= -ic4(F + l)i^i. (15) 



We solve Eq. (15) for t < and match the solution at t = 0. Note that 
the instanton solution F rapidly grows at t — with increasing A while it 
vanishes as i — > — oo. By using Eqs. (13) to (15), 

-^ln(F + l) + cc^(F + l)2^ = 0, (16) 

we find the solution, 

|ln(F + l) = c,{F + ir, (17) 

/ / \ l/2c 

F{t) = -1 + ^ + tanh(c^t - a2/2) . (18) 

The constants ai and a2 are to be determined using the boundary conditions 
F(0) = Fq and F(— oo) = 0. To find the value of F at t = we integrate Eq. 
(13) over (— e, e) and use Eq. (14) to obtain, 

Fi(0) ^ C3A. (19) 
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We then use Eq. (12) at i = to obtain 

FiOy^F^^-iX'^. (20) 
Using Eq. (20) and the condition that F{—oo) = 0, we find the values of 



a, = 1(1 ± ^1 - 4ci)2 and a, = -2tanh-^(Hi^^^^). 

We now find an approximate value of the PDF in the limit of A ^ oo 
which is expressed by a path-integral, using the saddle point method. First 
we compute the A-dependence of the action Sx and then proceed to determine 
the PDF tails. In the limit of A ^ oo, 5";^ becomes 



J dtc, (|fi(^ ln(F + 1) + C2(F + 1)^) j 

1 r _ 

-C1C4 J dtF^ 



+ -C1C4 

+ iX J dtciC3ln{F + l)S{t) 
= M^F'(0)+iAc3ln(F(0) + l) 

= -ikiX + ikiXln{-ik2X), (21) 
where the coefficient ki and k2 are defined as 

= (22) 

c 

k2 = (23) 

2C2 ^ ^ 

The tail of the PDF is then found by performing the A-integral in Eq. (3) 
by the saddle point method in the limit R ^ 00. It is later shown that this 

corresponds to A — cxd, 

P(R) ~ / ^^^iXR+ikiX-iki\\n{-ik2X) ^ \ £-^(^0) ^24) 

^ g-|^e«/'=i+i./(2.0_ (25) 
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Here, we have evaluated the A-integral using the extreme point /'(Aq) = 
and Ao = i-^e^^''^ of /(A) = —i\R + iki\ + ik2\\n\ by using the saddle-point 
method. Note that the first term (e '^z ) in Eq. (25) is the zeroth or- 
der term and (e'^/^^'^^^) is the correction term. The resulting PDF tail in Eq. 
(25) satisfies a Gumbel distribution [2lj- [25]. The Gumbel distribution is one 
type of an extreme value distribution that is the limit distribution of max- 
ima/minima of independent and identically distributed random variables. 
Thus the Gumbel distribution is one model of approximate maxima/minima 
of finite sequences of random variables. This is due to the fact that Taylor 
expansion of the exponential contains an infinite sum of all powers. 

We now show explicitly that systems with exponential non-linearities may 
have enhanced PDF tails compared to Gaussian distributions. In figure 1, a 
comparison of the PDF tails of fiux with an exponential non-linearity found 
in Eq. 25 with ki = 2.0 and k2 = 1.0 (black line), ki = 1.0 and k2 = 1.0 
(green line), ki = 0.5 and ^2 = 1.0 (blue line) and a Gaussian distribution 
(red line) are shown. We have normalized all the distributions as one-sided 
PDFs. The figure shows that the PDF tails with exponential non-linearity 
may be enhanced over a Gaussian distribution depending on the parameters 
ki and k2. It is important to note that as R oo the Gaussian distribution 
takes larger values than the Gumbel distribution the probability of taking 
such a large value of R may become extremely small (~ 10~^°). However, if 
the exponential non-linear interaction is weak and can be expanded in a trun- 
cated Taylor series up to second order the plasma sheath region is governed 
by a power-law non-linearity, a slightly different PDF tail ~ exp{— c(-R/-Ro)*} 
with s = {n + l)/m can be expected [22]. Here, n = 2 and m are the order 
of the highest dominant nonlinear interaction term and moments for which 
the PDFs are computed, respectively. 
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Figure 1: The PDF tails of flux in a system with exponential non-linearity 
with ki = 2.0 and k2 = 1.0 (black line), ki = 1.0 and k2 = 1.0 (green line), 
ki = 0.5 and k2 = 1.0 (blue line) and a Gaussian distribution (red line) are 
shown. 

Note that the usefulness of the found PDF tails is limited by the spatial 
structure of the specific solutions to the non-linear equations. The spatial 
structure determines the magnitude and sign of the coefficients ki and k2- 

Furthermore, note that to evaluate the accuracy of the Gumbel distri- 
bution in representing the PDF would require the exact computation of our 
path integral and A-integral in Eq. (3), which is a rather formidable task 
beyond the scope of the present paper. We can however obtain rough es- 
timate on errors in our predictions by considering the two major, leading 
order contributions that were neglected - the first one in using the saddle- 
point method for the path integral [Eq. (5)], and the second in the A-integral 
in Eq. (24). We note that we computed the PDFs keeping zeroth, first and 
second order terms in the A-integral and up to first order in using the saddle- 
point path integral. For instance, for R = 0.6 an error of the same order as 
the expansion in inverse aspect ratio ~ 30% (at i? = 1 the error is only a few 



percent) is found from the sub-leading order terms. Therefore, the Gumbel 
distribution is a reasonable estimate for i? > 0.6 — 0.7. The exact calculation 
of corrections to our model is however very involved, and will be addressed 
in future publications. 

In summary, in the present letter we have presented a statistical theory of 
turbulence and intermittency due to coherent structures in the sheath region 
in fusion devices. The use of coherent structures is motivated by various ex- 
perimental results in the sheath region that bursty events cause a significant 
transport which are linked to coherent structures. To investigate intermit- 
tency caused by these coherent structures, we have computed the PDF tails 
of 1st moment flux (e.g. density, potential) using the non-perturbative tech- 
nique called the instanton method. The PDF tails are found to be in the 
form of modified Gumbel distributions. Interestingly, the Gumbel PDF tails 
may be enhanced in comparison with a Gaussian distribution. It is impor- 
tant to note that PDF tails of flux were shown to be far from Gaussian 
although Gaussian forcing was used. This result could guide in interpreting 
experimental results. Specifically, the PDF tails derived here and in previous 
papers may directly be compared with experimental results using log-log or 
log-lin plots. Furthermore this model could elucidate the very complex in- 
termittent behavior at the edge of magnetic confinement systems. 
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